Destabilization of the Mg-H system through elastic constraints 
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Abstract 

We tune the thermodynamics of hydrogen absorption in Mg by means of elastic clamping. The loading isotherms measured by hydrogenog- 
raphy show that Mg films covered with Mg-alloy-forming elements, such as Pd and Ni, have hydrogen plateau pressures more than two orders 
of magnitude higher than bulk Mg at the same temperature. An elastic model allows us to interpret the Mg thickness dependence of the hy- 
drogen plateau pressure. Our results provide the basis for the development of new hydrogen storage materials with excellent thermodynamic 
properties. 
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Hydrogen is an attractive energy carrier for a future sus- 
tainable energy system. Compact hydrogen storage is how- 
ever still a scientific and technological challenge [1]. Stor- 
age of molecular hydrogen, both gaseous and liquid, requires 
high pressures or very low temperatures and is thus not en- 
ergy efficient. An alternative is to store atomic hydrogen in 
metal or complex hydrides. Storage in metal-hydrides allows, 
in principle, to achieve high volumetric and gravimetric den- 
sities and to reversibly operate at room temperature and at- 
mospheric pressure. However there are thermodynamic and 
kinetic limitations associated with the chemical reactions in- 
volved in hydrogen absorption and desorption. 

An optimal metal-hydride system has a hydrogen equilib- 
rium pressure of 1 bar at room temperature. Due to weight 
constraints the interest of the scientific community has turned 
to lightweight hydride-forming elements such as Li, B, Na, 
Mg and Al. Mg can store up to 7.6 mass% of hydrogen, form- 
ing the ionic MgH2 compound, but shows poor kinetics of 
hydrogen absorption/desorption and a hydrogen equilibrium 
pressure of 1 bar at the far too high temperature of 573 K. The 
absorption and desorption kinetics can be enhanced by reduc- 
ing the size of the Mg grains via ball milling |2] hence short- 
ening the hydrogen diffusion length, or by addition of proper 
catalysts pT. Many efforts are currently dedicated to mod- 
ifying the thermodynamics of the Mg-H system in order to 
promote hydrogen dissociation at lower temperatures. Some 
destabilization has been observed upon reduction of particle 
size [4] and alloying of Mg with other elements IH 0], al- 
though this generally implies a reduction in the total hydrogen 
storage capacity. Theoretical works suggest the possibility to 
reduce the hydrogen desorption temperature in both MgH2 
nanoparticles with grain size smaller than ^^1.3 nm |7p and 
in MgH2 layers with thicknesses below few unit cells [8]. 

In this letter we show that the thermodynamic stability of 
the Mg-H system can be drastically modified via elastic con- 
straints. This is a direct consequence of the long range H-H 
interaction in metals [l^. In a metal with free surfaces 
two hydrogen atoms feel an attractive interaction, while in a 
perfectly clamped sample the interaction becomes repulsive. 



Clamping can thus be used to tune the equilibrium hydrogen 
pressure to the levels required for specific applications. To 
prove the validity of our approach we prepared two series of 
capped Mg thin films: 1) to explore the role of chemical bind- 
ing we deposited Mg thin films capped with different transi- 
tion metals: Ni, Pd, Ti, Nb and V; 2) in order to study the 
thickness dependence of the elastic effect we prepared a se- 
ries of Pd-capped Mg films of different thicknesses. The ef- 
fect of clamping is strongly dependent on the chemical nature 
of the transition metal used as a cap-layer and will be inter- 
preted with a simple elastic model. By choosing the appro- 
priate sample geometry we are able to increase the hydrogen 
plateau pressure more than 200 times with respect to bulk Mg. 
Hydrogen absorption in our samples is measured by means of 
optical spectroscopy and Hydrogenography a novel tech- 
nique that allows to measure Pressure-Optical Transmission- 
Isotherms (PTI's). In a PTI the amount of light transmitted 
by a thin film is measured as a function of increasing pres- 
sure at constant temperature: when the metallic Mg films load 
with hydrogen the Mg-MgH2 metal-insulator transition leads 
to an abrupt increase in the amount of transmitted light. From 
the Beer-Lambert law, ln(T/TM) oc end, the logarithm of 
the optical transmission, T, normalized for the transmission 
of the film in its metallic state, Tm, is proportional to the hy- 
drogen concentration in the film, ch and to the film thickness, 
d. Measuring PTI's is therefore in all respects equivalent to 
measure standard Pressure-Composition-Isotherms [11]. The 
films used in this study are deposited in a UHV chamber 
(p < 10^^ Pa) by RF/DC magnetron sputtering on glass 
and quartz substrates. All films are covered with Pd to pre- 
vent oxidation and catalyze hydrogen absorption. Optical re- 
flection spectra are measured in a PerkinElmer Lambda 900 
diffraction grating spectrometer, with energy range between 
0.5 and 6.5 eV, during (de)hydrogenation at room tempera- 
ture and H2 pressures up to 10^ Pa. PTI's are recorded in 
an optical cell that allows to vary the hydrogen pressure from 
10^^ to 10^ Pa at constant temperatures, between 313 and 
573 K. Further details on the hydrogenography experimen- 
tal setup can be found in Ref. [6]. In Fig. [T]the PTI of 



1 



a 2x[Ti(10nm)Mg(20nm)] film covered with 20 nm of Pd is 
shown. In the pressure and temperature ranges explored Pd 
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FIG. 1: (color online) Loading PTI measured at 333 K for 
a 2x[Ti(10nm)Mg(20nm)]Pd(10nm) sample deposited on glass. 
In the inset a cartoon of the sample geometry: glass- 
substrate/Ti/Mg/Ti/Mg/Pd. 



and Ni do not load with hydrogen. Ti, Nb and V layers on the 
contrary load at the beginning of the pressure ramp but, due 
to the metallic nature of their hydrides, produce very small 
optical effects compared to Mg. The observed plateaus are 
therefore only due to the formation of MgH2 upon hydro- 
gen absorption. The isotherm shown in Fig. [T] exhibits two 
clear and distinct plateaus of similar width. Since the width 
of the plateau is proportional to the thickness of the material 
the two plateaus must originate from the two distinct Mg lay- 
ers. Furthermore, to obtain such a well defined double-plateau 
we have to assume that the two Mg layers are elastically dis- 
connected [12] and have different site energies for hydrogen 
absorption. To prove these hypotheses we measure the opti- 
cal reflection of the 2x[Ti(10nm)Mg(20nm)]Pd(20nm) sam- 
ple deposited on quartz, while slowly increasing the hydrogen 
pressure. Optical spectra are measured through the transpar- 
ent quartz substrate at near-normal incidence. In Fig. [2^ the 
measured spectra are shown for three different steps of load- 
ing: 1) as-deposited metallic sample, 2) at intermediate load- 
ing and 3) fully loaded sample. The calculated spectra shown 
in Fig. I2J3 are generated with SCOUT |13] by taking into 
account the volume expansion occurring upon hydrogenation 
of Mg and Ti and including a rough interface between the 
quartz substrate and the first Ti layer. In the simulation the 
intermediate step is due to the complete loading of one of 
the two Mg layers while the other remains metallic. Com- 
parison with the experimental spectra shows unambiguously 
that the Mg layer sandwiched between two Ti layers absorbs 
hydrogen at a lower pressure than the one capped with Pd. 
This is rather counterintuitive since the "Ti-sandwiched Mg" 
lies on the bottom of the sample and hydrogen has to dif- 
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FIG. 2: Measured (a) and simulated (b) optical reflection measured 
through the quartz substrate for a 2x[Ti(10nm)Mg(20nm)] film cov- 
ered with 20 nm of Pd at different stages of loading. In the simula- 
tions two possible states are taken into account for the intermediate 
loading: the "loading from top" and "loading from bottom". The di- 
electric functions of Ti, Mg, TiH2, MgH2 and Pd used in the simula- 
tions are taken from Ref. [IJ], [IJ], yj], [17|] and [3] respectively. 



fuse through the top Mg layer, which remains in the diluted 
ct phase, in order to reach the bottom one. For this mech- 
anism to take place the two Mg layers have to differ with 
respect to their thermodynamic properties and they must be 
elastically disconnected by the presence of a Ti layer in be- 
tween. This remarkable behavior is the result of the follow- 
ing H loading sequence. Already at very low H2 pressures Ti 
forms TiH2. The consequent 25% volume expansion, together 
with the positive enthalpy of solution of Mg and Ti, leads to 
disconnected interfaces between Mg and TiH2 layers. The 
Ti-sandwiched Mg layer is then effectively elastically discon- 
nected from the surrounding. However, since Pd alloys with 
Mg and loads at higher pressures, this reasoning is not appli- 
cable to the Ti-Mg-Pd block. When the top Mg layer expands 
upon H absorption it feels the Pd elastic constraint, result- 
ing in a higher hydride formation enthalpy and consequently 
in a higher plateau pressure. To demonstrate the validity of 
this interpretation we considered two series of samples: a) 
Mg thin films capped with different transition metals, and b) 
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Pd-capped Mg thin film with different thicknesses. Fig. [3] 
shows the PTI's for 5 different samples with the following ge- 
ometry: Ti(10nm)Mg(20nm)X(10nm)Pd(10nm) with X = Ni, 
Pd, Ti, Nb and V. Due to the presence of a Ti layer between 
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FIG. 3: (color online) Effect of cap-layer: PTI's measured at 333 
K for Ti(10nm)Mg(20nm)X(10nm)Pd(10nm) samples deposited on 
glass with X = Ni, Pd, Ti, Nb and V. The dashed line is the pressure at 
which coexistence of a and /3 phases begins to appear upon hydrogen 
absorption in bulk Mg | 18.1 . In the inset: schematic visualization of 
the elastic model used to interpret the data. 



Mg and the substrate we can assume that the only clamping 
effect is due to the top capping layer Clearly two behaviors 
are distinguishable: 1) Mg alloy forming elements such as Ni 
and Pd have a strong effect on the thermodynamic properties 
of Mg, leading to a much higher plateau pressure than the one 
measured for Mg films of equal thickness capped with Ti, Nb 
and V; 2) the latter elements, which are immiscible with Mg, 
effectively behave like "scissors" and lead to "quasi-free" Mg 
layers with elastically disconnected interfaces. Assuming that 
the origin of the destabilization observed is only elastic in na- 
ture, we can build a simple model consisting of two layers - 
a Mg and a cap-layer - glued to each other as shown in the 
inset in Fig. [3] Hydrogen atoms absorbed by the Mg layer 
isotropically expand the lattice of the metallic host, while the 
Mg feels a compressive stress in the x and y directions due to 
the cap layer In the z direction there is no stress built up by 
the clamping but we have to take into account the strain along 
the z direction due to the Poisson's ratio, v. The two lay- 
ers are in vacuum and no effect of the substrate is taken into 
account, basically assuming that the bottom Ti layer is a per- 
fect scissor By taking into account only elastic constraints we 
can calculate the effective volume expansion of the Mg layer 
upon hydrogen absorption. Furthermore, the enthalpy of hy- 
dride formation has a simple volume dependence through the 
bulk modulus, B HI: dA7?/dlnV^ = -BMgVij, where Vh 
is the partial molar volume of hydrogen in Mg. We can then 
derive an expression for the plateau pressure of the capped 



magnesium layer, Mg*, with respect to free bulk magnesium, 
Mg°: 
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where E — \RTVyijV^, V^ig is the molar volume of Mg 
and E = 3B{1 — 2v) is the Young's modulus. The equation 
can also be rewritten as [ln(p*/p'')]^^ = / + Sdyi^ with 
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According to the model the equilibrium pressure for the 
capped Mg has a straightforward dependence on the thickness 
of the Mg layer, dug- Fig. |4] shows the PTI's for 5 samples 
with different Mg thicknesses: Ti(10nm)Mg(znm)Pd(40nm), 
with z = 10, 15, 20, 30 and 40 nm. Two effects are clearly 
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FIG. 4: (color online) Effect of Mg thickness: PTI's measured at 333 
K for Ti(10nm)Mg(znm)Pd(40nm) samples deposited on glass with 
z = 10, 15, 20, 30 and 40 nm. In the inset: Mg thickness dependence 
of the equilibrium pressure. The dashed line is a linear fit to the 
points. 



visible: a) the width of the plateau, w, is proportional to dMg 
as expected from the Beer-Lambert law, b) the plateau pres- 
sure decreases with increasing Mg thickness as predicted by 
the model. A 10 nm thick Mg film has an equilibrium hy- 
drogen pressure of 2.8-10^ Pa at 333 K, more than 200 times 
higher than bulk Mg at the same temperature (P333 — 12 
Pa ill). The plateau width is measured by dividing each 
isotherm into three regions: the solid solution (a phase) be- 
fore the plateau, the coexistence region (a + f3 phases) and 
the fully hydrogenated region (/? phase) after the plateau and 
fitting each part with a straight line: the width is then given 
by the distance, in units of ln(r/TM), between the intersec- 
tion points. With the same geometrical construction we obtain 
the equilibrium pressure, as the pressure in the middle of the 
plateau. Ideally one would expect ui — > for ^Mg 0, how- 
ever, although the plateau width is found to be proportional 
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to the Mg thickness, a linear fit shows a relation of the type 
w = A{dMg — do), with do — 6 nm. The presence of an inter- 
cept indicates that the total amount of Mg undergoing a phase 
transition to MgH2 is thiner than expected and that we can 
assume for all samples that ~6 nm of Mg are "lost" due to al- 
loying with the top Pd layer In the comparison with the model 
we will therefore take into account an effective Mg thickness 
d = c?Mg — df). The formation of an interfacial Mg-Pd alloy 
is directly responsible for the strong clamping effect observed 
and is consistent with the difference in width observed in the 
two plateaus of Fig. [T] corresponding exactly to 6 nm. The 
inset in Fig. |4]shows the thickness dependence of the plateau 
pressure of Pd-capped Mg layers: plotting [ln(p*/p'^)]^^ ver- 
sus the effective Mg thickness we find a linear dependence as 
predicted by the model, with intercept / ~ 0.16 and slope 
S = 7.6 • 10^"^ nm^^. Substituting into equations |2] the lit- 
erature values for the structural parameters of Mg, Pd and H 
(Vh ^ 2.24 cm^/mol, ^Mg = 13.97 cm^/mol, vug = 0.29, 
j/Pd = 0.39, Eug = 45.2' GPa, Ep^ = 126 GPa) and tak- 
ing T = 333 K and dpd = 40 nm, we obtain the theoret- 
ical values Ith = 0.27 and Sth = 2.1 • 10"^ nm"^ The 
agreement between experiment and theory is surprising when 
considering the strong approximations made in the model: no 
substrate is taken into account, the only interaction is elastic, 
hence neglecting surface energy contributions, and the layers 
are assumed to behave as perfect elastic bodies. The rather 
large S value obtained from the fit can be explained consid- 
ering the alloying effect taking place at the Mg/Pd interface 
and assuming that only a fraction of ~10 out of 40 nm of 
Pd has a significant clamping effect on the Mg underneath. 
The discrepancy in the / values could reflect the fact that we 
consider the Young's modulus of metallic Mg: taking into ac- 
count MgH2, which is stiffer than Mg, would lead to a closer 
agreement between experiment and theory. The elastic model 
developed for a double-layer can be easily extended to 3D Mg 
nanoparticles embedded in a skin of hard material. Assuming 
Mg particles of 10 nm of diameter and a skin of 2 nm thick- 
ness, the model predicts an increase in hydrogen equilibrium 
pressure of more than 2000 times with respect to bulk Mg at 
the same temperature. The effect of clamping described in the 
present work is not only important in the perspective of hydro- 
gen storage applications, but it is crucial in understanding the 
thermodynamics of all hydrogen-absorbing capped thin films, 
from switchable mirrors [20.1 to electrode materials for batter- 
ies 112 ill , hydrogen detectors ll22ll multilayers and superlattices 
Q. 

In conclusion we have shown the possibility to tailor the 
thermodynamics of a metal-hydrogen system by means of 
elastic constraints. Thin films of Mg are used as a model 
metal-hydrogen system. Elements immiscible with Mg be- 
have like scissors, while Mg-alloy-forming elements exert a 
clamping effect that leads to huge increases in hydrogen equi- 
librium pressures. Clamping arises as a consequence of al- 
loying between Mg and the capping layer and its effects can 
be understood on the basis of a simple elastic model. The 
possibility to tune the thermodynamics of a metal-hydrogen 



system by elastic means offers attractive new possibilities for 
compact hydrogen storage. 
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